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Abstract 

Using a recently proposed factorization hypothesis for semi-inclusive hard 
processes in QCD, one can study, in principle, the diffractive production of 
the Standard Model Higgs boson at LHC using only, as input, ep diffractive 
hard-processes data of the type recently collected and analyzed by the HI and 
ZEUS collaborations at HERA. While waiting for a more precise and complete 
set of data, we combine here the existing data with a simple Pomeron-exchange 
picture and find a large spread in the Higgs boson production cross section, 
depending on the input parametrization of the Pomeron's parton content. In 
particular, if the Pomeron gluon density f g /p(P) is peaked at large j3 for small 
scales, single diffractive events will represent a sizeable fraction of all produced 
Higgs bosons with an expected better-than-average signal-to-background ra- 
tio. 
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Searching for the Higgs boson is one of the main goals of future hadron col- 
liders. The main search strategy for this elusive particle rests on the produc- 
tion via the dominant gg — > H fusion process, and on the observation of the 
subsequent decays into ZZ and ZZ*, for the mass ranges m H > 2m z and 
130 GeV < van < 2m z, respectively (the experimental signature being four 
charged leptons from the Z, Z* decays), or of the rare decay into two pho- 
tons, H — > 77, for the intermediate mass range 90 GeV < ran < 130 GeV. The 
search for the Higgs boson in the intermediate mass range is believed to be 
very difficult |TJ. The process pp^gg + X^H + X^^yy + X suffers indeed 
from the large irreducible background pp — > 77 + X and from the reducible 
background pp — > jet + jet + X, jet + 7 + X, where the jets fake photons |@||. 

Additional information on the parton-level initial state (i.e. distinguishing gg, 
qg or qq induced processes) in a selected class of Higgs boson production events 
may help to reduce the background. For instance, the irreducible two-photon 
background is mainly induced by a quark-ant iquark pair and is therefore 
suppressed if one is able to bias the sample so as to favour initial state gluons 
(see the brief discussion of backgrounds at the end of the paper). In hard 
diffractive processes the proton stays essentially intact and its valence quarks 
go straight into the leading final proton. What recoils against the proton- 
proton system is naturally expected to be gluon-rich, at least in comparison 
to unbiased (i.e. diffractive plus non-diffractive) events. 

Diffractive hard processes are usually described in the framework of QCD by 
introducing an effective flux for the exchanged object with vacuum quantum 
numbers (for short, called the "Pomeron" (P) in the sequel) and by parametriz- 
ing the parton densities of the Pomeron itself M. There is some experimen- 
tal evidence in support of the theoretical expectation that the Pomeron is a 
gluon-rich object ||||]7]]. Tagging leading protons with a large momentum ||, 
or observing gaps in rapidity as an experimental sign for diffractive processes, 
may allow to exploit additional information about the hard scattering process 
and to enrich the sample of gluon-initiated events. 

In this paper we study the production of the Standard Model Higgs boson 
in hard diffractive processes at LHC, using experimental fits to the Pomeron 
parton densities from groups working at HERA. We point out, however, that, 
for a study of hard diffractive events, this (partly model dependent) theoret- 
ical framework is not really needed. A direct parametrization of the gluon 
content of a proton fragmenting into a leading proton would be sufficient and 
is desirable in order to facilitate forthcoming studies of hard diffraction. The 
corresponding theoretical framework of fracture functions has been developed 
recently PJTOl. The fracture function Mp h (x, z, Q 2 ) gives the joint probability 
distribution for an observed hadron h (e.g. a leading proton) with a specific 
momentum fraction z in the proton fragmentation region, and a parton % (e.g. 
a gluon) of momentum fraction x = (3(1 — z) initiating the hard process. The 
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functions ([3, Q 2 , x-p) introduced in [11 ,12] are closely related to fracture 
functions in the same way that F 2 (x, Q 2 ) is related to the usual parton densi- 
ties. We wish to encourage an experimental analysis in terms of such a direct 
parametrization and thus independent of the "Pomeron exchange picture" 
which, being non-perturbative, is beyond present theoretical control. 





(a) 



(b) 



(c) 



Fig. 1. Generic diagrams for the total Higgs boson production cross section (a), 
the single diffractive case (b) and the double diffractive case (c). The quark in the 
triangle loop is the top quark. Single lines stand for the incident proton and leading 
outgoing proton, and double lines for the fragments of the proton and Pomeron. 



Figure 0a shows the standard gg fusion graph for Higgs production via a 
top quark loop, giving the total cross section for pp — > gg + X — > H + 
X in leading order, the process being diffractive or not. Either one or both 
protons may stay essentially intact, giving rise to single (Fig. |I|b) and double 
(Fig. |l]c) diffractive processes. For a soft gluon content of the Pomeron, the 
double diffractive process has been studied in |L3|], and the formalism of Regge 



theory has been applied in |L4] . We base our study on the framework for hard 
diffraction developed in [15|, also employed in |13|], and do not make use of 
Regge theory, but instead assume the picture of the Pomeron as an exchanged 
object carrying vacuum quantum numbers, with a parton content that can be 
measured and parametrized. 



The flux /p/ p (xp, t) of Pomerons in the proton is given by the parametrization 
of Donnachie and Landshoff |16,17|]. Actually, since the rapidity-gap criterion 
does not distinguish a leading proton from some other diffractively excited 
i n should be increased relative to the Donnachie-Landshoff value. 



state, 



Equivalently, one can follow the practice of experimental groups to stick to 



the flux of Refs. 16.17 but allow the overall normalization of the Pomeron's 



parton densities to exceed the bounds imposed by a naive momentum sum 
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rule. For the parametrizations of the Pomeron's gluon density^ we use results 
from fits by the HI ||[7j and ZEUS f|l2] collaborations (parametrizations 
Hla, Hlb, Hlc and ZEUS0) and moreover employ recent results by Gehrmann 
and Stirling |H| (parametrizations GS1 and GS2). In addition, to facilitate a 
comparison with ||13| , we also include two parametrizations with a soft gluon 
content (parametrizations S and SE, without and with scale evolution, re- 
spectively). The parameters of the input distributions are shown in Tab. [TJ. In 
the cases where the parametrizations were not available at arbitrary scales, 
the evolution has been carried out with the standard leading-order Altarelli- 
Parisi equations with heavy flavour thresholds at the single quark masses^ 
and Aqcd = 200 MeV for 4 flavours. 



Figure [| shows the various gluon densities and the quark singlet distribution of 
the Pomeron at a scale of 100 GeV. Two of the HI parametrizations (Hla and 
Hlc) form a lower and upper bound of all gluon parametrizations considered 
here; only the ZEUS parametrization is smaller at small values of (3. As we 
will see later on, the Pomeron is probed mainly at large /3, and therefore 
the large spread of the parametrizations of about one order of magnitude 
translates into a corresponding variation in the calculated cross sections. The 
quark singlet distribution is better constrained by the measurement, 
and so the distributions have a smaller spread at large /3, where the present 
measurements are most sensitive. 



1 We do not include the model of Buchmiiller and Hebecker [[rj] in our study, 
although it gives a first principles QCD formulation of diffractive events without 
employing the concept of the Pomeron, because it is not clear how to define the 
gluon density for diffractive events in this case. The concepts of "rotation in colour 
space" and "colour-compensating soft gluon exchange" of this model must probably 
be applied to the virtual top quark in the triangle loop. Because of two incident 
gluons, there should be an additional constraint related to the colour matching. 

2 The HI parametrizations are based on fits of F® , whereas the ZEUS fits also 
take into account constraints from the photoproduction of jets. In the latter case, 
if only the information from F^ is used, it is possible to have a consistent set of 
parameters with a large gluon content as well. To be definite, we use the central 
values given in ||. 

3 The evolution for the fit of the parameters in the F% analysis of the HI col- 
laboration has been done with Nf = 3 flavours, and by adding charm via the 
photon-gluon fusion process. In our case, the evolution span is sufficiently large in 
order to justify the treatment of the charm and bottom quarks as massless flavours. 
The influence of changing Nf from the values used in our procedure to a fixed num- 
ber of 3 flavours is to increase the gluon density, because this decreases the number 
of quark flavours into which the gluon can split. As a result, the cross section for 
the single diffractive case increases by about 10-20%, and the one for the double 
diffractive case by about 20-40%. 
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Table 1 

Parametrizations of the Pomeron's parton content. The input distributions are given 
by Pfg/wW) = AP B (1 - (3) c for the gluon and 0f q/P (0) = Df3 E {l - (3) F for N light 
quark flavours q, q at the scale ^Lq. (3 is the momentum fraction of the gluon in the 
Pomeron. We have also included the line types that are used in the plots later on. 
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Fig. 2. Parametrizations of the Pomeron's gluon density (a) and singlet quark dis- 
tribution (b). The line types are given in Tab. @. 

The cross section for Higgs boson production via gg fusion [^(J is given by 



0" = j y 91(^,^)92 (jf^^J THa ° ( 



H 



III 



top , 



The quantity th is denned by m 2 H / Eq M , where Ecu is the centre-of-mass 
energy of the collider. An explicit expression for the parton-level cross section 
cr o{ m H / m top) can be found, e.g., in |2T|. For the total cross section and in 



the single diffractive case, both or one of the gi{£,,fi 2 ) are the gluon densities 
fg/p{£,, A* 2 ) of the proton, respectively, which we choose to be the GRV leading- 
order parametrization f22fl . The other g^s, in the single and double diffractive 
case, are given by the convolution of the Pomeron flux factor fp/ p (xp,t) and 
the parton densities of the Pomeron f g /p({3,n 2 ): 

0t(f,/^ 2 ) = / di / — - fp/ P (xp,t) fg/v — ,/i 2 . (2) 

J J Xp \ Xp I 

— oo £ v ' 



The upper limit of the xp-integration is given by the condition that the 
Pomeron's momentum fraction of the proton does not exceed a certain frac- 
tion xp max , in order not to invalidate the parametrization used for the Pomeron 
flux factor. Moreover, there is a kinematical limit depending on the momen- 
tum transfer t = (p' — p) 2 , where p and p' are the momenta of the incident 
and outgoing proton, respectively, so that 



7(f>*) = max { f> min 



1m 2 p 



(3) 



m p being the proton mass. In this expression we have neglected terms of the 
order of m p /E CM and \/— t/-BcM> reflecting the ambiguity in the massive case 
(m p ^ 0, t ytz 0) of the definition of the frame- dependent momentum fraction 
variable xp. 

For the numerical evaluation we have chosen Eqm = 10 TeV, m top = 180 GeV 
and xp max = 0.1. Increasing the centre-of-mass energy to 14 TeV results in an 
increase of all cross sections by about a factor of 2. The factorization scale 
[L and the renormalization scale are set to the Higgs boson mass. We use 
the leading-order running strong coupling constant with Aqcd = 200 MeV 
for 4 flavours. The dependence of the cross section on the Higgs boson mass 
w,h is shown in Fig. |^a for the single diffractive caseQ and in Fig. |]b for 
the double diffractive case. The shape of the mass dependence for the sin- 
gle diffractive case is similar to the non- diffractive case, although the cross 
section drops faster for increasing Higgs boson mass. In the mass range of 
90 GeV < rriH < 130 GeV, the single diffractive cross section for the Hlc 
parametrization is about 25% of the total cross section, whereas most of the 

4 We have subtracted the double diffractive cross section from the single diffrac- 
tive one, so what is actually plotted is the contribution where the proton which is 
modelled by the GRV gluon density does not stay intact. 
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Fig. 3. Dependence of the single diffractive (a) and double diffractive (b) Higgs boson 
production cross section on the Higgs boson mass. The full line is the total cross 
section pp — > gg+X — > H+X. The other line types related to various gluon densities 
of the Pomeron are defined in Tab. [J. 



other parametrizations (excluding Hla) studied here give a fraction of the or- 
der of 10%. Except for the Hlc parametrization, the double diffractive cross 
section is very small, less than about 0.1 pb. Varying the gluon content of the 
ZEUS parametrization within the bounds given in || yields a variation of 
about ±60%. A comparison of the cross sections for the parametrizations S 
and SE shows that the evolution of the parton densities gives, in the double 
diffractive large effect; however, this obviously depends on the input 

scale chosen. It should be noted that the diffractive cross section is not yet 
well constrained by the measurements of the Pomeron structure function. For 
ran = 100 GeV, the spread is about a factor of 14 in the single diffractive case 
and about a factor of 600 (!) in the double diffractive case. 

Differential distributions for the momentum fraction variables x-p, (3, £pr = 
xp(l — P), £jv.D and £d = x-p (3 are shown in Fig |j. Here £p r is the momen- 
tum fraction relative to the proton momentum carried by the remnant of the 
Pomeron, ^^d is the momentum fraction of the gluon from the fragmenting 
proton in the single diffractive case and £p is the momentum fraction of the 
gluon from the Pomeron, relative to the momentum of the incoming proton. 
The momentum fraction f3 of the gluon in the Pomeron is large, and the 
distribution flda/dfl peaks at about 0.2. The distribution of the momentum 
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x dff/dx x dcr/dx 




Fig. 4. Distributions xda/dx for various variables x in the single diffractive (a) 
and double diffractive (b) case for a Higgs boson mass of 100 GeV and the Hlb 
parametrization of the gluon density. The integrated cross section a is normalized 

to 1. The variable x stands for xp [ ], f3 [ ], xp(l — (3) [■■■], £nd [ ]> 

£,D [—■—]■ The ZEUS parametrization gives similar results. 



fraction xp of the Pomeron in the proton, being cut off by the explicit upper 
limit xp max , has large contributions at large scp, although the Pomeron flux 
behaves like 1/xp. The reason for this behaviour is that the invariant mass 
squared of the produced system in the hard scattering process is fixed and 
quite large, i.e. th = m 2 H /EQ M = 10~ 4 for mjj = 100 GeV and Ecm = 10 TeV. 
The momentum fraction £ to be supplied by each of the protons is therefore 
of the order of 10 -2 (cf. Fig. f|). Since the gluon density of the Pomeron is de- 
creasing rapidly for increasing f3, the Pomeron itself must carry a comparably 
large fraction of the proton momentum. This fact has the consequence that 
the momentum fraction £p r carried by the "Pomeron remnant" relative to the 
proton momentum is, on the average, much larger than 10 -2 . The rapidity 
gap which would be observed in principle is between the beam direction and 
the Pomeron remnant. Assuming naively that the latter is simply a jet with 
longitudinal momentum £p r p and transverse momentum of the order of the 
transverse momentum of the scattered proton, multiplied by (1 — /3), then it 
turns out that the distribution of the pseudorapidity r\ of this jet peaks at val- 
ues of the order of r\ = 7. This means that hardly any gap will be observable 
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in real experiments^, even if fragmentation effects are taken into accountQ. 
Instead, for diffractive events, under the assumption that the scattered proton 
escapes detection, the signature will be a comparably small hadronic activ- 
ity in the forward direction, with the unobserved scattered proton carrying a 
large momentum fraction^. In the double diffractive case, the measurement 
of the invariant mass of the two-Pomeron system is in principle possible by 
measuring the energy of the two scattered protons in forward spectrometers 
|| . If the Pomerons fragment, then the measurement of the invariant mass of 
the two-Pomeron system probably gives no constraint on the invariant mass 
of the two-gluon system, because the gluons carry only a small fraction of the 
Pomeron momentum. There may, however, be a statistical correlation of xp 
and £d. If the Pomerons do not fragment, then a measurement of the momenta 
of the scattered proton gives the mass of the produced Higgs boson. The cross 
section of this case, the so-called elastic part, has been evaluated in [|TJ] in the 
framework of Regge theory and is found to be of the order of 0.1 pb, which is 
of the same order of magnitude as what we get for the total double diffractive 
case, including possible fragmentation of the Pomerons. This situation should 
be clarified. 

In order to make more reliable predictions for hard diffractive scattering at 
pp colliders, the gluon content of the Pomeron has to be determined much 
more precisely. The present experimental analyses cover a range from about 
Q = 3 GeV to Q = 10 GeV, whereas the factorization scales to be employed 
for pp collider physics reach up to 1 TeV. The gluon density at large scales is 
not very well constrained by the analysis of at small scales, because it 

contributes in 0(a s ) only, and because of the large evolution span. Additional 
constraints by studying heavy quark and jet production at HERA and possibly 
at the Tevatron are needed to improve the situation. 

The viability of the present study rests on the assumption that the simple 
factorization picture of eqs. (|I|) and (0) applies. The analysis of a toy model 
shows that the concept of factorization may break down for certain diffractive 
processes [ffijl , and thus the expression of the cross section as a convolution 
of parton densities with a mass-factorized parton-level scattering cross section 
ceases to be valid in principle, although it may still be a good approximation 
in reality. This issue may be studied at present high-energy colliders. The pro- 

5 It is to be expected that the same argument applies also to the production of 
other systems with large invariant mass, such as jets and heavy quarks. 

6 The Pomeron remnants will be dragged towards the hard process because of the 
missing colour connection to the incident proton. However, this should not really 
influence the discussion, give or take one unit of rapidity. 

7 If the mass of the Higgs boson is assumed to be very small, of the order of 10 GeV, 
then it turns out that the xp and (3 distributions are much broader, and that the 
pseudorapidity distribution of the Pomeron remnant jet extends to smaller values 
of r\. 
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duction of heavy quarks or jets at large pr in diffractive events in pp collisions 
is closely related to Higgs boson production studied in this paper, in the sense 
that the produced system has a large invariant mass and thus the momen- 
tum fractions probed in the proton and Pomeron are large. A comparison of 
the corresponding cross sections from the Tevatron with theoretical predic- 
tions would give a hint of what can be expected in the case of Higgs boson 
production at the LHC. 



Finally, we wish to make some comments on the background to the process 
studied here. As said before, the idea is to exploit the fact that the Pomeron 
is a gluon-rich object, and thus background processes that are quark initiated 
should be suppressed. We briefly consider two cases. One is the decay of the 
Higgs boson into t + t~, proposed in [24j]. It has been shown in |25[] that there 



is an overwhelming background from tt production that renders this channel 
unfeasible. Unfortunately, tt production proceeds mainly via gg fusion, and 
so this background is expected to be large in the diffractive well. It 

might be the case, however, that the tt production cross section turns out 
to be smaller, because of the large invariant mass of the tt pair. The other 
case we wish to comment on is the decay H — > 77. As mentioned in the 
introduction, there are two types of backgrounds, the irreducible one from the 
production of photons via qq, gg fusion and bremsstrahlung, and the reducible 
one, where jets fake photons. The latter background and the bremsstrahlung 



contribution can be reduced considerably by photon isolation cuts p6 |. The 
direct production of photon pairs in qq and gg fusion has been studied in 
PTfl , where it has been shown that the qq process dominates. This background 
might thus be reduced in the case of diffractive events. It seems to be worth 
while to amend these qualitative remarks about the background by a thorough 
quantitative studyQ. 

To summarize: we have studied the production of the Standard Model Higgs 
boson at the LHC in single and double diffractive processes. Depending on the 
parametrization of the gluon content of the Pomeron, the production cross 
section for the single diffractive case may be sizeable and reach up to 25% 
of the total Higgs boson production cross section, if the Pomeron is gluon- 
dominated at large f3 for small scales. However, the present parametrizations 
for the gluon density of the Pomeron are not yet sufficiently precise to allow for 
a reliable prediction. Due to the large mass of the Higgs boson, the diffractive 
events studied in this paper have only a "small" rapidity gap between the 
leading proton and the Pomeron fragments, the latter having typically a very 
small transverse, but an extremely large longitudinal momentum. Besides a 
leading proton detected in a forward detector, the experimental signature will 



We would like to remark that an effective gluon density for hard diffractive pro- 
cesses may easily be implemented in existing programs by parametrizing the con- 
volution integral in eq. (|2|), thus replacing a parton density by a fracture function. 
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be a comparably small hadronic activity in the forward direction. 

We wish to thank J. Philips for discussions about the parametrizations of the 
Pomeron's parton densities, L. Lonnblad for clarifying comments on rapidity 
gaps in electron-proton scattering, J. Terron for a cross check of the evolu- 
tion of the parton densities, L. Trentadue for conversations about fracture 
functions, and Z. Kunszt for discussions about the physics of the Higgs boson. 
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